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Abstract: Water availability, which enables plant growth and animal activity, regulates dryland ecosystem 
function. In hyper-arid ecosystems, rain cannot support vascular plant growth. Therefore, hyper-arid 
vegetation is restricted to the lower topography, where runoff accumulates. Typically, food resources 
originating from areas of dense vegetation are dispersed across the desert floor, enabling animal life in 
areas lacking vascular plant growth. However, certain regions, such as the hyper-arid upper topography, 
may be devoid of plant-derived food resources. The present study examined arthropod activity in the 
upper topography of a hyper-arid desert, in comparison with arthropod activity in the lower topography. 
Pitfall traps were utilized to compare arthropod activity along unvegetated ridges with activity in parallel, 
vegetated riverbeds. Surprisingly, the study revealed dense arthropod communities in the barren upper 
topography. Arthropods collected in the upper topography represented 26% of total arthropod 
abundance. In addition, the overlap between arthropod identity in the ridges and wadis (i.e., riverbeds) was 
low, and certain arthropods were strongly affiliated with the ridges. The upper topographic communities 
included high numbers of silverfish (Zygentoma: Lepismatidae), malachite beetles (Pséloderes), and 
predatory mites (Acari: Anystidae), and these arthropods were present at various life stages. It remains 
unclear how arthropod communities can persist in the unvegetated upper topography of the hyper-arid 
study area. These results raise the possibility that other food sources, independent from vascular plants, 
may play a significant role in the life history of hyper-arid arthropods. 
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1 Introduction 


Drylands are defined as regions with an aridity index of less than 0.65; the aridity index is 
the ratio of mean annual precipitation to potential evapotranspiration (Barrow, 1992). Drylands 
cover 40% of the Earth's land area and support over two billion people (White and Nackoney, 
2003). These areas face the double threat of climate change and land degradation (Giorgi et al., 
2001; DI Giovanni and Palombini, 2002; IPCC, 2007; Trenberth, 2011). Drylands include dry 
subhumid, semi-arid, arid and hyper-arid lands. Hyper-arid lands are defined by an aridity index 
below 0.05 (Barrow, 1992). This quantitative threshold reflects distinctive shifts in the spatial 
distribution of vegetation and associated trophic systems. Unlike other drylands, hyper-arid 
lands feature large 
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areas beyond the dry limit of plant life, where minimal precipitation and high rates of evaporation 
prevent vascular plant growth (Noy-Meir, 1973). While extensive research has focused on arid 
ecosystems, hyper-arid ecosystems are understudied. 

In dryland ecosystems, water availability regulates plant growth and drives primary 
productivity. Water availability in dryland ecosystems is regulated by runoff redistribution, or the 
flow of water from "source" areas to "sink" areas (Noy-Meir, 1973). On the community scale, 
impermeable soils act as water "sources", while permeable soils act as water "sinks" (Pueyo et al., 
2013). On the landscape scale, water is redistributed in different land morphology; sloped or 
elevated areas act as water "sources" and level or low-lying areas act as water "sinks" (Cammeraat, 
2002; Mayor et al., 2008). These patterns of water availability regulate the vegetation abundance 
in dryland ecosystems (Schwinning et al., 2004). Much like dryland vegetation, dryland animal 
abundance also depends on water availability. Many animals in dryland have developed 
independence from dietary water (Ward, 2009). Nonetheless, these animals rely on food resources 
that originate from primary producers. Accordingly, the fauna appearance in dryland often reflects 
the distribution of vegetation and water availability (Parker and Pianka, 1975; Meserve and Glanz, 
1978; Dean et al., 1999). Such spatial correlations between primary and secondary productivity 
are found in many terrestrial ecosystems (Valavanis et al., 2004; Holtgrieve et al., 2009; Bregović 
and Zagmajster, 2016). 

In hyper-arid ecosystems, water availability is all the more critical to plant and animal 
distribution. Hyper-arid landscapes feature large areas too dry to support plant life. As a result, 
hyper-arid plant communities are restricted to isolated pockets on the desert floor, or to the lower 
topography (Monod, 1954; Zohary, 1962; Kassas and Girgis, 1965). The concentration of 
vegetation in the lower topography (namely dry riverbeds and canyons, termed "wadis"), reflects 
the central role of land morphology and runoff in the distribution of water availability. Within the 
lower topography, soil properties play a central role in hyper-arid plant distribution (Abd 
El-Ghani and Amer, 2003; El-Bana and Al-Mathnani, 2009). The ridges of hyper-arid ecosystems 
are completely bare from vascular plants since moisture in the ridges soil is not sufficient to 
support them (Zohary, 1962; Noy-Meir, 1973; Fossati et al., 1999). 

As in other dryland ecosystems, the fauna of hyper-arid ecosystems is often distributed in areas 
where vegetation grows, either as a result of trophic dependency, or because of shared preferences 
for cooler and less dry environments (Rebelo and Brito, 2007; Nawaz et al., 2011; Hackett et al., 
2013; Alatawi et al., 2020). However, certain arthropods, such as silverfish (Zygentoma) and 
darkling beetles (Tenebrionidae), are not restricted to vegetated areas, and instead are widely 
distributed across the desert floor. Even so, these arthropods still rely on nutritious organic 
detritus, produced near vegetation and dispersed by wind or animal migration (Louw and 
Hamilton, 1972; Watson, 1989). Some hyper-arid ecosystems, such as the core of the Atacama 
Desert, Chile, lack both plant life and any form of plant-derived food resources (Ruhm et al., 
2020). The only organisms known to survive such harsh environment are specialized microbial 
communities (Warren-Rhodes et al., 2006). Conditions are similar in the upper topography of 
many hyper-arid ecosystems, where water is too scarce to support vascular plant growth. Due to 
the lack of available food resources, hyper-arid upper topographic areas have never been targeted 
by ecological studies. 

The present study examined the spatial distribution of arthropods in the hyper-arid ecosystems 
in southern Israel. Arthropod activity in the steep upper topographic areas, lacking vascular plant 
life, was compared with that in the lower topographic areas, where runoff accumulation supports 
plant life. Arthropod activity was used as an indicator of ecological activity, because arthropods 
occupy various ecological niches and respond to fine-scale spatial variation in environmental 
conditions (Kremen et al., 1993; Liu et al., 2016). Vegetation and topography were used as 
indicators of environmental conditions. The researchers did not record soil moisture, because rain 
in the study area is too scarce to create reliable differences in soil moisture. The researchers 
hypothesized that arthropods would be found only in the vegetated range of the study area. 
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2 Materials and methods 


2.1 Study area 


The study was conducted in the highlands near to the west of the Arava Institute for 
Environmental Studies (AIES), which belongs to the Kibbutz Ketura of the Arava Valley, 
southern Israel. The Arava Valley is a desert basin that extends along the border of Jordan and 
southern Israel, from the Dead Sea to the Red Sea. The climate of the southern Arava Valley is 
hyper-arid, with an annual mean temperature of 30°C, mean annual precipitation between 30 and 
50 mm, and mean annual evaporation of 3000 mm (resulting in an aridity index of 0.010—0.017). 
During the last century, temperatures have increased while precipitation has decreased (Ginat et 
al., 2011). The valley is bounded by the Edom Mountains to the east and by highlands to the west. 
The elevation of the valley ranges from 400 m below sea level to 210 m above sea level 
(Goldreich and Karni, 2001). The study was conducted in three sites each with two landscape 
units: one ridge and one riverbed adjacent to the ridge (Fig. 1). The ridges were bare of any 
vegetation, while the riverbeds were consistently vegetated with both shrubs and acacia trees. 


2.2 Sampling methods 


Four linear plots were established at each site, half in the wadis and half on the ridges, for a total of 
12 plots. Each ridge plot was paired with a riverbed plot to create 6 blocks, 2 at each site (Fig. 1). 
Distance, change in elevation and slope degree are provided in Table 1, along with photos of the 
terrain (Fig. 2). Each ridge plot was established on a separate peak, at least 100 m away from 
vegetation. Each riverbed plot was established in the center of the riverbed, within 5 m of vegetation, 
and at the point closest to the corresponding ridge plot. At each of the 12 plots, 4 dry pitfall traps 
were placed at 1-m intervals. Cups with a diameter of 7 cm were used. In the riverbed, the traps 
were established on loose sand and gravel, while in the ridges, a pickaxe was required to break the 
rocky soil and establish the traps. Over a period of 6 weeks, 48 traps were sampled 4 times. 
Arthropods that could not be identified in the field were brought to the lab. Other arthropods were 
photographed on-site and identified to the highest recognizable taxonomic unit (RTU) (Ward and 
Stanley, 2004; Bird et al., 2020). 
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Fig. 1 Location and topographic features of the study sites in the Arava Institute for Environmental Studies 
(AIES), southern Israel 
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Table 1 Site conditions at each block, including distance between the ridges and wadis, change in elevation 
between the ridges and wadis, and slope degree in the ridges and wadis 


Block Distance (m) Change in elevation (m) Wadis slope (°) Ridges slope (°) 
A 70.7136 10.42416 17.08256 16.66176 
B 84.4296 10.24128 5.65228 6.21687 
C 55.7784 7.31520 22.94079 10.18023 
D 95.4024 8.71728 15.08426 15.94099 
E 137.7696 18.40992 8.15679 15.46925 
F 173.7360 24.59736 15.75063 4.91531 


Fig.2 Photos showing terrain in the ridges (a) and wadis (b) of the plots 


2.3 Data analysis 


Arthropod activity was calculated by dividing the number of individuals trapped by the number of 
days that the trap was open (Rubinstein et al., 2013). Arthropod activity, Simpson diversity and 
species richness in the ridges and wadis were compared by paired-sample t-test, with each subsite 
as a replicate. Dissimilarity between the ridges and wadis was calculated using the Whittaker and 
Bray-Curtis indices, and evaluated against zero (no difference between the ridges and wadis) by 
one-sample f-test, with each subsite as a replicate. These analyses were completed in R version 
4.0.0. 

Ants, beetles and soil mesofauna respond to specific ecological factors, necessitating separate 
t-tests. Ants respond to the pheromones of other ants. Beetles in the Arava Valley are 
morphologically diverse, but the vast majority is wandering detritivores. Soil mesofauna includes 
arthropods between 2 mm and 100 um in diameter (Swift et al., 1979). These arthropods are 
highly abundant and may be ecologically linked to microorganisms. Without justification for 
further subdivision of the data, all other arthropods were grouped together for t-tests. Flying 
insects, including Lepidoptera, Diptera, Aphididae and winged Hymenoptera, were excluded from 
the analysis. Pitfall traps are not suitable for flying insects, therefore, flying insects found in the 
traps were not considered. 

Multivariate analysis linked RTUs and locations using CANOCO 5.0. Principle component 
analysis (PCA) was used to cluster locations by similarity between their community compositions 
(ter Braak and Smilauer, 2015). Environmental variables were tested for clusters of locations 
using redundancy detrended analysis (RDA). The RDA ordination tests if the locations (ridges 
and wadis) of the arthropod assemblies are different from each other by taking location 
(topography) as an environmental variable. Interactive-forward selection was completed including 
only significant variables step by step. 
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3 Results 


3.1 Activity, richness, diversity and dissimilarity 


The majority of arthropods (74%) were trapped in the wadis. Arthropod activity was significantly 
greater in the wadis than in the ridges for ants and soil mesofauna, but not for beetles or other 
arthropods (Table 2). Simpson diversity was significantly greater in the wadis than in the ridges, 
for all groups except ants (Fig. 3; Table 2). Species richness was significantly greater in the wadis 
than in the ridges, for all groups. The Whittaker and Bray-Curtis dissimilarities were significantly 
greater than zero for all groups. 


1.0) (a) CORidge E Wadi 
os ; 
3 0.64 E 
5 a 
o 
5 0.4] 3 
a A 
n 0.2 kr 
0.04 : : - : 
S Site 1 Site 2 Site 3 Site 1 Site 2 Site 3 
oO 
2, 405 (c) 0.9) (d) MM Bray-Curtis CI Whittaker 
B 35 0.84 
& 304 0.7] See KAX] 
3 2064 S : SeSe 
a 25: 5 REX ` eases 
g 2 g 25] RS neces 
a 70 £04) R bS =- SeS 
B15 3 LS] bS? =: SSS] 
=e = 0.34 <> PO. : L] 
£ = SSI eels Kx) 
= 104 0.24 So RI B Seq 
“ oS) RX) ` Xs 
2 ISS Peg SSS 
5 5] 0.1) SS Pee: KS] 
Š o 00 KZ 2020x ERE [SS 
< Site 1 Site 2 Site 3 Site 1 Site 2 Site 3 


Fig. 3 Simpson diversity (a), species richness (b), arthropod activity (c) in the ridges and wadis, and Whittaker 
and Bray-Curtis dissimilarities (d) 


Table 2 Results of t-tests (df=5) activity, Simpson diversity, species richness and Whitaker and Bray-Curtis 
dissimilarities for ants, soil mesofauna, beetles and other arthropods 


t-test Ants Soil mesofauna Beetles Other arthropods 
Activity (abundance/days P=0.0104 P=0.0118 P=0.2178 P=0.6796 
elapsed) 1=3.3294 13.2155 t=1.4093 t= -0.4381 
P=0.2792 P=0.0111 P=0.0213 P=0.0092 
Simpson diversity 
1=0.2134 1=3.2674 t=2.7044 1=3.4383 
P=0.0021 P=0.0215 P=0.0014 P=0.0014 
Species richness 
t=5.0000 t=2.6968 t=5.4515 t=5.5000 
P=0.0033 P=0.0184 P<0.0001 P<0.0001 
Whittaker dissimilarity 
t=4.4779 t=2.8274 t=11.3830 t=11.2330 
P<0.0001 P=0.0003 P<0.0001 P<0.0001 
Bray-Curtis dissimilarity 
t=21.8110 t=7.6950 t=34.4910 t=10.6280 


3.2 Multivariate analysis 


The results of PCA show a clear distinction between the communities of the ridges and wadis (Fig. 
4). RDA shows that several RTUs are associated with the ridges (arrows to the right), but most 
RTUs are affiliated with the wadis (Fig. 5). The constraint analysis of RDA showed that 
topography was significant and explained 36.7% of the variance (n=12, pseudo-F=5.8, P=0.004), 
whilst sites were not significant (P=0.350; Fig. 5). 
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Fig. 4 Principle component analysis (PCA) plot of the recognizable taxonomic unit (RTU). Filled circles 
indicate ridges and open circles indicate wadis. Numbers by each circle indicate the number of species in each 
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Fig. 5 Redundancy detrended analysis (RDA) plot showing affiliation of recognizable taxonomic units (RUTs) 
with either ridges or wadis. Filled circles indicate ridges and open circles indicate wadis. 
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3.3 Community composition in the ridges and wadis 


Numerous soil mesofauna were captured in the pitfall traps, including elongate bodied springtails 
(Collembola: Entomobryomorpha) and mites of the Anystidae family. The springtails were mostly 
restricted to the wadis, but the mites were distributed evenly across the ridges and wadis (Fig. 6a). 
Predators were scarce in the ridges, with the exception of jumping spiders (Salticidae; Fig. 6b). 
Ants were also scarce in the ridges, with the exception of Cataglyphis nigra, which was 
repeatedly captured at one subsite (Fig. 6c). Detritivores were numerous in both the ridges and 
wadis (Fig. 6d). The detritivore population of the ridges was dominated by silverfish (Zygentoma: 
Lepismatidae) and malachite beetles (Psiloderes). A few herbivores and fungivores were collected 
in the wadis, but none were collected in the ridges (n=7). 

Certain arthropods were more numerous in the ridges than in the wadis (Fig. 6). These included 
the jumping spiders, silverfish and malachite beetles, as well as one species of ant, and one 
species of darkling beetle. In addition, various life stages of certain arthropods were observed in 
the ridges, namely silverfish and rove beetles (Staphylinidae). In fact, 16 of the 18 rove beetle 
larvae were collected in the ridges. These larvae are either predators or detritivores. 
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Fig. 6 Community composition in the ridges and wadis. (a), soil mesofauna; (b), predators; (c), ants; (d), 
detritivores. A-F are 6 blocks. 


4 Discussion 


The results of this study clearly reveal two striking conclusions. First, although vegetation is 
absent from the ridges, life in the ridges is abundant, making up 26% of total arthropod 
abundance (Fig. 4). Second, communities in the ridges are predominantly distinct from those of 
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the riverbeds below, and not merely a reduced abundance of riverbed communities. This is 
evident from the clusters revealed by PCA, the affiliation of certain arthropods with the ridges, 
and the low degree of overlap between the identity and activity of arthropods in the ridges and 
wadis (Figs. 2 and 3). Certain arthropods, namely silverfish (Zygentoma: Lepismatidae), 
malachite beetles (Psiloderes), and predatory mites (Anystidae), were consistently observed at 
high densities, across all sampling sites in the ridges (Fig. 4). It remains unclear how these 
arthropods complete their life cycle within the harsh conditions of the ridges. 

Traditionally, hyper-arid arthropods are thought to depend on vascular plants. Herbivores are 
thought to be restricted to vegetated areas, and the survival of arthropods beyond the limit of plant 
life is explained by wind-blown organic detritus, produced near vegetation (Louw and Hamilton, 
1972; Crawford and Taylor, 1984; Watson, 1989). It is quite possible that the results of the present 
study are merely one more example of this allochthonous detritus paradigm. However, growing 
evidence has linked hyper-arid arthropods to lithophilic autotrophs, such as lichens (Shachak et 
al., 1987; Lalley et al., 2006). In the past, these links have only been demonstrated in areas of 
dense lichen coverage, and in proximity to vascular plants. However, lithophilic autotrophs such 
as cyanobacteria are well known from the driest environment of the Earth, and may be ubiquitous 
on sufficiently stable hyper-arid substrates (Warren-Rhodes et al., 2006; Pointing and Belnap, 
2012). Previous studies of allochthonous detritus in deserts are largely restricted to sand dunes, 
where the unstable substrate prevents the development of robust microbial communities 
(Crawford and Taylor, 1984). In contrast, the rocky cliffs of the study area provide an ideal 
environment for the development of cryptic autotrophs such as lichen and cyanobacteria. Many of 
the arthropods observed in ridges, such as silverfish and darkling beetles (Tenebrionidae) are best 
known as detritivores, but these groups are also known to feed on lichens (Furniss and Carolin 
1977; Wessels et al., 1979). 

If the arthropod communities of the ridges do in fact rely on a local food source, this would 
bring into question the life history of many hyper-arid arthropods, which are thought to be 
dependent on plant-derived food resources. Such a discovery would provide a novel perspective 
on the functioning of hyper-arid ecosystems, and might significantly impact other areas of 
dryland ecology. Currently, however, the ecology of the ridges remains unclear. Due to the 
unexpected results of the study, additional research will be required to shed light on the ecology 
of the ridges, including surveys of lithophilic autotrophs. Surveys of ridges farther removed from 
riverbeds would clarify the relationship (or lack thereof) between the arthropods of the ridges and 
the vegetation of the wadis. Genetic barcoding of silverfish gut contents could further clarify the 
food sources available to arthropods of the ridges. 


5 Conclusions 


The unvegetated, hyper-arid ridges of the study area support novel arthropod communities, which 
differ from the communities of neighboring vegetated areas. The dominant arthropods in the 
ridges, including silverfish (Zygentoma), malachite beetles (Psiloderes), and darkling beetles 
(Tenebrionidae), are best known as detritivores. Arthropods were represented at various life stages 
in the ridges, and predators such as jumping spiders (Salticidae) and certain mites (Anystidae) 
were also well represented. It remains unclear what food sources are available in the ridges, 
where conditions are too harsh to support plant life. However, the abundance of arthropods in the 
ridges raises the possibility that local food sources, independent from vascular plant life, could 
play a key role in ridges communities. Further investigation of trophic systems in the harsh 
environments of the ridges may provide novel perspectives on the life history of hyper-arid 
arthropods, and the structure of hyper-arid ecosystems. 
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